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We study two similar spin ladder systems with the ferromagnetic leg coupling. First model includes 
two sorts of spins, s = 1/2 and s = 1, and the second model comprises only s = 1/2 legs coupled by 
a "triangular" rung exchange. We show that the antiferromagnetic (AF) rung coupling destroys the 
long-range order and eventually makes the systems equivalent to the AF s = 1/2 Heisenberg chain. 
We study the crossover from the weak to strong coupling regime by different methods, particularly 
by comparing the results of the spin-wave theory and the bosonization approach. We argue that 
the crossover regime is characterized by the gapless spectrum and non-universal critical exponents 
are different from those in XXZ model. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Gb 



I. INTRODUCTION 

The strongly correlated systems in one spatial dimen- 
sion (ID) attracted an enormous theoretical and experi- 
mental interest last decade. The ID fermionic and spin 
systems were recognized long ago as the useful theoretical 
models, where the interaction effects are very important 
and at the same time are subject to rigorous analysis. "'^^^ 
The experimental discovery of the systems of predomi- 
nantly ID character inspired the renewed interest to this 
class of problems. Among the experimental realizations 
of the ID systems one can mention the Bechgaards salts, 
carbon nanotubes, copper oxides spin ladders and purely 
organic spin chain compounds.^ 

While the physics of purely onc-dimensional objects, 
or chains, is well understood,-^'^ the spin ladders are 
still under intense investigation. Even the isolated spin 
chains reveal a variety of unusual phenomena, including 
the Haldane gap, spin-Peierls transition, magnetization 
plateaus. The ladders, consisting of a few coupled spin 
chains are generally much richer in their behavior, and 
pose additional theoretical problems. 

The interest to the problem of ladders may be traced 
back to the earlier attempts to construct the continous 
representation of spin 5 = 1 variable in ID out of s = 1/2 
quantities.i^- The methods elaborated in these studies are 
now widely used in the analysis of the ladder systems. 

The basic model for the spin chain, the antiferromag- 
netic (AF) Heisenberg s = 1/2 chain, is thoroughly stud- 
ied by various methods.^ This may be one of the rea- 
sons, why a majority of the theoretical papers discussing 
the spin ladders are now confined to the treatment of 
quantum spin s = 1/2 with antiferromagnetic interac- 
tion along the legs. Spin ladders with different spins or 
with a ferromagnetic leg exchange attracted much less 
attention. 

The spin chains and ladders consisting of different 
spinsi and with the AF leg exchange were considered re- 
cently in^iSiLgiS. Particularly, the ferrimagnetic chains 



with alternating spin-1/2 and spin-1 were discussed 
there. It was shown that, contrary to the case of equal 
spins, the uncompensated spin value in a unit cell leads 
to the gap in the spectrum and the appearance of the 
long range magnetic order ii£ 

The spin-1 /2 ladder with ferromagnetic exchange along 
the legs was studied iniiiS. A rather rich phase diagram 
was demonstrated, depending on the details of the mag- 
netic anisotropics for the leg and rung couplings. 

In the present paper we study the mixed spin (5=1, 
s — 1/2) spin ladder with the ferromagnetic exchange 
along the legs and antiferromagnetic interaction on the 
rungs. In the absence of the rung coupling, the individual 
chains show the ferromagnetic long-range order (LRO), 
and the ground state is classical. The spectrum and the 
ground state energy is well described in terms of the lin- 
ear spin-wave theory.^ We show that the inclusion of 
the antiferromagnetic rung coupling drives the system 
into the quantum regime, understood in terms of the 
AF Heisenberg spin-1/2 chain. In this regime the mag- 
netic LRO is absent, and the spatial correlations show 
the power-law decay. Note that the uncompensated spin 
in a unit cell leads to the absence of the gap in low-lying 
excitation spectra. 

This primary observation is interesting on its own, be- 
cause two limiting cases of the model allow the asymp- 
totically exact solutions with gapless spectra. Hence our 
further motivation to study the crossover region starting 
from both limiting points. 

We discuss the regimes of weak and strong rung cou- 
pling for different types of the leg exchange anisotropy. 
The consideration is somewhat complicated by the ab- 
sence of the established routines for our case. The 
bosonization, an extremely useful tool in dealing with 
AF s=l/2 systems, does not fully work here on two rea- 
sons. First is the existence of two sorts of spins in a unit 
cell, and another is the ferromagnetic leg exchange. 

Hence we supplement our study by the consideration of 
a similar model, written entirely for s = 1/2 but with the 
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FIG. 1: A ladder of spins S = 1, s = 1/2 



modified rung couplings. Using these models and com- 
paring diS^erent approaches, the spin-wave theory and 
bosonization, we arrive at the unified description of the 
ferrimagnetic spin ladders. Particuarly, we discuss the 
spectrum and the correlations and observe the crossover 
from the weak to strong coupling regime. An attention 
is paid to a subtler point in the bosonization procedure, 
a seemingly unstable Gaussian effective action near the 
ferromagnetic point. 

The rest of the paper is organized as follows. We dis- 
cuss the mixed spin model in strong and weak rung cou- 
pling regime in Sec. ^1 The spin s — 1/2 ladder with 
"triangular" rung exchange is introduced and analyzed 
in Sec. IIIII The existence of different order parameters 
in a system is discussed in Sec. IIVI The discussion and 
conclusions are in Sec. Ivl 



doublet, described by a spin s = 1/2 variable a2i- We 
show below, that the effective interaction between these 
CTi is antiferromagnetic for the above choice J > and 
hence the situation is mapped onto a well-known problem 
of s = 1/2 Heisenberg antiferromagnet. One docs not 
expect the long-range order at D = and the spectrum 
is linear ~ Jk. The correlation functions for this case 
are described below. 

The intermediate situation, J ~ Jj^, _D 7^ 0, is harder 
to analyze and we present several ways to discuss it be- 
low. 



A. Strong rung coupling 

Consider first the case of strong perpendicular cou- 
pling, J± ::§> J. Taking first J ~ D ^ 0, one sees 
that each pair of spins s, S" is characterized by a total 
spin j = s -I- S and by the rung energy Ej = J±sS = 
1-^-1 + 1) -s(s + l)-S'(S'+l)). The lowest state here 
is doublet, j = 1/2, the first excited state is quadruplet 
j — 3/2, with £^3/2 — Ei/2 = 3J±/2. The wave function 
of a multiplet is 
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(2) 



II. 



MIXED SPIN LADDER 



We investigate the properties of a ladder system, con- 
sisting of two sorts of spins, s = 1/2 and S' = 1, arranged 
in a checkerboard manner. The unit cell comprises four 
spins and the Hamiltonian is 



i 

i + Si,2i+lS2,2i+l) (1) 

i 

with the first subscript (e.g. 1 in Si^2i) labeling the leg, 
and the second one denoting the site on it (odd or 
even). Below we mostly consider the case of the AF 
rung coupling, Jj^ > 0. The overall ferromagnetic ex- 
change along the chains allows the uniaxial anisotropy, 



J; 



jy = J > 0, J||= 



J + D > QAD\ <^ J. In what 



follows, we consider also a useful generalization to higher 
spins, s ^ 1 with the difference S — s = 1/2 kept fixed. 
The whole consideration is done for zero temperature. 

First let us briefly describe two simple limiting cases. 
At -D = 0-t- and Jj_ = 0, we have two ferromagnetic 
chains, which possess the long-range magnetic order. The 
spectrum is quadratic at small wave vectors, w ~ Jq^. 
We put the lattice spacing a to unity everywhere except 
the Sec. ImDl 

In the opposite limiting case, ^ J, the ground 
state of two spins on the rung, say, Si^2i and 82,22 is 



with Clebsch-Gordan coefficients C-'™ ■ The inter- 
action along the chain is considered now as a perturba- 
tion. The formula Q and consideration below are ap- 
plicable for larger spins as well. In this more general 
case of s = S — 1/2 > 1/2, the operators s", 5"" act 
within the lowest doublet and connect the doublet with 
the quadruplet, but the direct transitions to the higher 
states, j = 1/2 ^ j > 5/2, are absent. One can check 
that the corresponding matrix elements are given by 



(j = l/2,m|s"|j-l/2,m') = 

(j = l/2,m|5"|j = l/2,m') = ^0„„ (3) 

with cr^„j/ the Pauli matrices. For s = 1/2 this reads 
1 



6 



(4) 



It shows that if we consider only the projection of the 
spin operators onto the lowest doublet, then the above 
Hamiltonian corresponds to the AF Heisenberg spin-1/2 
model of the form : 

^e// = ^ Jeff^i<^i+l . (5) 



where 



J^ff{Ji_ -> 00) = 2J|f s(5 + l)/9. (6) 
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Below we will refer to this estimate of J"f 
For s = 1/2 we have J"ff (oo) = 2J^/9. 

Let us next consider the role of the higher states on a 
rung. We do not list here the matrix elements of the spin 
operators S", s" for the transitions j = 1/2 ^ j = 3/2. 
We only note that they are proportional to \/ s{S + 1) 
and are the same for S" and s", except for the sign. 
The second-order correction in J|| between the adjacent 
rungs, labeled by 1 and 2 below, can be written as 



2Ei/2 — 2i?3/2 



(7) 



with factor 4 coming from the above property of coinci- 
dence of the matrix elements. Noting that 



s{S+l) 
9 



we find that the second order in J|| results in © with 
the renormalized value of effective interaction : 



j'^f^{oo)J2ff{co) 



(9) 



with J"jj from For the above parameters J" and 
s = 1/2 one has 



J, 



2 SJiJ + D) 
9 729Jj_ 

r ' 729 Jj_' 



(10) 
(11) 



which particularly means that the relative value of 
anisotropy decreases with decreasing Jj_. The perturba- 
tion theory is expected to break down when the correc- 
tion in is comparable to the first term in (|5J|, which 
happens roughly at ^ "^if ■'^^te that this crite- 
rion also corresponds to the point where the bandwidth 
induced by Je// becomes comparable to the separation 
~ Jl between the quadrupet and doublet. At larger 
the low-energy dynamics is described by the AF spin- 
one-half XXZ model (jSJ which is exactly solvable. 

In the isotropic case, I? = 0, the spectrum is linear, 
Lo = f^e//!?! and the correlations are of the form 



i-\-ni 



+ (-l)"|n|- 



with the omitted minor logarithmic corrections 

Taking into account the matrix elements we find 
that the leading asymptotes in the isotropic, 5 = 1/2,5" = 
1 case are : 



\2n\'' 
4|2n-f- 
16|2n| 



(12) 



Thus the long-range ferromagnetic order is absent, but 
the correlations are ferromagnetic and slowly decaying. 
In addition, there is a subleading sign-reversal asymp- 
tote and modulation depending on the spin value (s — 
1/2,1). The correlations between the spins in different 
chains are slowly decaying antiferromagnetic ones, e.g. 
{^1 i+2n+i) ~ —\2n + 1|~^. The above form of the 
correlations is obviously unchanged as long as J± < J. 



B. Weak rung coupling, spin- wave analysis 

Let us consider next the opposite limiting case, when 
the leg exchange dominates and J± can be considered 
as perturbation. First we explore the spin-wave formal- 
ism for the easy- axis anisotropy, D > 0. It was shown 
recently^ that in case of ferrimagnetic system in ID the 
spin-wave description gives very good estimate for the 
ground state energy and on-site magnetization. Perform- 
ing the standard Dyson-Maleyev expansion. 
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(13) 



we write for the magnon Green function in the linear 
spin-wave theory (LSWT) approximation 

$^ = (4,5j',c_fc,d_fc) 
G,,{k,t) = -te{t){[^,{k,t),^Uk,0)]) (14) 



/s(2J| + Jj_)-w VsSJ-fk Vs5J_i 

VJSJjk S{2J^^ + J±_) - uj VJSJjL 

V7SJi_ si2JI + Jj_) + uj VsSJjk 

\ VsSJjL VsSJ-fk S{2^ + Ji_)+ujJ 



(15) 
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Here 7fc = 2 cos k. The spectrum consists of doubly de- 
generate acoustic and optical modes. The optical mode 
has the energy ^ ^/sSJ for all wave vectors and the 
acoustic branch at small k is 



2sS 
s + S 



^J{2D + Jk^){2D + Jk^ + 2Jj_) (16) 



so instead of the quadratic spectrum of purely FM case, 
we have an approximately linear spectrum at small en- 
ergies Ek S sJ±. The contribution of the zero-point fluc- 
tuations into the average on-site magnetization can be 
estimated for Z) ^ Jj^ <C J as follows 



S - [Sli 



^/JjJH^JljD). 



(17) 



It means that the spin-wave approximation fails when 
the latter quantity is of order of s. It happens roughly 
at J_L > J {hi J / D)~'^ , which corresponds to the 
crossover point 7* inii. Apart from the logarithmic fac- 
tor, this estimate agrees with the above value ^ '^^'^11 ' 
obtained in the large limit. For lower the system 
shows the long-range order. 

It might be instructive to consider the FM rung cou- 
pling J± < 0. In this case all the branches of the spec- 
trum are gapful, the lowest modes are 



£l,fc 



£2,k 



2sS 
s + S 

2sS 
s + S 



{2D + Jk^), 

{2D + 2\J^\ + Jk^), (18) 



Note that at the isotropic point D = 0, the long-range 
FM order and hence the applicability of the LSWT is lost 
for any J± > 0. The spectrum in this case is gap- 
less and linear with the spin- wave velocity ~ s\/ JJ_l, the 
limit of s ~ 5 ^ 1 is assumed. We know that increas- 
ing J_L we should eventually recover the effective model 
characterized by spinon velocity ^ s^J. These two 
velocities match again at J±^ ^ s^J. 

Concluding this section, we also present the LSWT 
results for the lowest branches of dispersion for the case 
of easy-plane anisotropy, D < 0. For the AF sign of the 
exchange J± one has 



£2,fc 



2sS 
s + S 
2sS 



y/j^H2\D\ + jW+2Jl), 
^{2\D\ + Jk^){2J^ + Jk^), (19) 



s + S 

while for the FM exchange Jj^ < we obtain 
2sS 



S2,k 



S + S 

2sS 
s + S 



y/jk^2\D\ + Jk^), 



(20) 



v/(2| J±| + Jk-^){2\D\ + 2| Jx| -I- Jfc2), 



In this case LSWT is formally inapplicable, but the above 
formulas might be useful for a comparison with further 
results. 



Summarizing, we observe that while the LSWT cannot 
treat the correlations correctly and formally is inapplica- 
ble in the absence of the LRO, it provides simple and 
reasonable formulas for the excitation spectra in a com- 
plicated one-dimensional ladder. We illustrate this point 
below by discussing the spin-1/2 ladder, where a rigorous 
description of the low-energy action is available. 



C. Equations of motion 

This subsection is devoted to the macroscopic equa- 
tions describing the behavior of spins in two coupled 
chains. To make our calculations more transparent, we 
define new operators on a rung as superpositions of two 
spin operators 



Q2i — ^2.2i + ^1.2i^ — '^'2,21 ~ ^l,2i; (21) 

'^2i+l — "-"1,21+1 T" *2,2'i+l; ^2i+l — '^l,2i+l *2,2i-|-li 

satisfying the following commutation relations: 

[R?,Q^] = ^^cp^Rl (22) 



with CajS-Y totally antisymmetric tensor ; the operators 
Q" and R" commute at different sites. 

The commutation relations H22|l define the 
SU{2)<^SUi2) = 50(4) group with two Casimir 
operators given by 



(23) 



The Hamiltonian H = + H± takes the form 



2 



5:(Q?-R^) = :^5:Q^ + c.t (24) 



The equation of motions for operators Q and R are given 
by 



(25) 



These equations correspond to the well-known Bloch 
equations for precession of the magnetic moment of fer- 
romagnets (antiferromagnets) . Taking into account that 
the operators (^") are different on two different 
sublattices corresponding to odd(even) sites, we use the 
properties 



fR 



2 yoal 



0, 



[Q^,i?f] = -[R^,i?r]=zea/3^{i?fQZ} 
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and adopt a symbolic notation 

ie^p^{R^Q2} - I (R. X Q, - Q, X R,)" . 

Then we obtain the foUowing system of coupled Bloch 
equations for the isotropic J|" 

dtQi = Y (Qj X Q,,±i - R, X Rj±i) 
9tR, = Y (Qj X R,,±i - Rj X Qj±i) + 

+ ^ (Q, X R, - R, X Q,) . (26) 
Taking the continuum limit here, one has: 

m = A^pQ X ^ - B^pR X J— J— 

( 9^R \ 

dtR = BapQ X — — A^pR X ^ + 

\ OXaOXp OXaOXp J 

+ C(QxR-RxQ). (27) 

where for asymmetric two-leg ladder a = (3 = x and 
Axx — Bxx — J\\0''^/'2, C = J±/'^- The lattice index i is 
omitted. Being re-written in terms of densities of mag- 
netic moments, equations H27|l corresDond to the gener- 
alized Landau-Lifshitz equations for the dynamic SO {4) 
group. Note that the operator Q represents a total spin 
on the rung, and the operator R does not have a simple 
interpretation. 

We introduce the densities of magnetic moment char- 
acterizing two sublattices 

O = |ER'' T = |E(-1)'R'' (28) 

with J\f being total number of rungs in the ladder. 

Obviously, M and N represent the uniform and stag- 
gered magnetizations along the ladder, whereas T and 
O can be interpreted as "staggered rung" magnetization 
in the uniform and staggered channels along the chain, 
respectively. The different ordered phases are charac- 
tererized by nontrivial values of M,N, 0,T or some 
combinations of them whereas in disordered phase these 
quantities are equal to zero. For example, the ordered 
FM phase is characterized by M 7^ 0. In Neel phase 
N ^ 0. In ferrimagnetic phase both M ^ and N ^ 0. 
Therefore, the information about the order parameter is 
necessary for deriving the macroscopic Landau-Lifshitz 
equations from the microscopic Bloch equations. 

We point out that the scalar product of two spins on 
a rung 

S1S2 = I (Q' - R') 
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FIG. 2: A triangular ladder of spins s = 1/2 



may also be considered as a local order parameter (see 
discussion in the Section III). 

Eqs. H26|l are the central result of this subsection. Upon 
the assumption of the certain order parameter (s) they 
lead to macroscopic "equations of motion" . Applying a 
standard routine^ one can show that these equations are 
equivalent to LSWT treatment and reproduce magnon 
dispersion laws discussed in the previous subsection. The 
detailed analysis of the Bloch equations for asymmetric 
50(4) ladders will be presented elsewhere. 



III. TRIANGULAR s = 1/2 LADDER 

One may regard spin 1 as a ground state triplet of two 
spins 1/2 coupled ferromagnetically. Our model assumes 
that this triplet is also FM coupled to other spins 1/2 on 
a chain. Therefore instead of the FM chain of n spins 
1 and n spins 1/2, one may consider 3n spins 1/2. The 
model we propose is 



i 

+ E (®2,3i + S2,3i+l) 

i 

+ (Sl,3i+1 + Sl,3i+2) S2,3i+2] (29) 

with the above choice of J|". The model (I29II is not 
equivalent to the previous one, eq. Indeed, the exact 
mapping of Q to (|29|l would include the strong trimer- 
ized isotropic FM in-chain exchange, at those links which 
form the bases of the triangles in Fig. [3 In that case one 
would first consider the triangles and then couple them 
to each other, fully restoring the consideration of the pre- 
vious Section. We show below that the model (|29|l with 
uniform value of the rung exchange has two advantages. 
First, it is equivalent to at J_l ^ oo and second, it 
is easier tractable in the opposite case ^ 0, since the 
exact form of the low-energy action is available for the 
uniform J,?. 
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A. Strong rung coupling 

We consider a case when the AF exchange J± is much 
larger than J|p first for the isotropic — this 
case a main block is a triangle formed by two rungs, and 
the coupling of triangles is a perturbation. The Hamil- 
tonian for the triangle is 



1 increase with 



Ha = Jjl{si + 82)83 - Jii (8182 - 1/4) 



(30) 



The structure of the energy levels is as follows. The 
term J|| groups the spins on one leg, 81 and 82, into a 
triplet |T) and a singlet \S). The singlet does not couple 
to 83 and results in a total doublet denoted as \D0) with 
the energy Ejjq = J||. The triplet |T) of zero energy cou- 
ples to 83 with the formation of doublet \D1) and quadru- 
plet \Q). The corresponding energies are i?_Di = —J± and 
Eq = J±/2. For J|| 00 the state \D0) is unimportant, 
and 81, 82 act as one spin S = 1. At the same time, the 
low-energy sector of the problem is associated with the 
doublet \D1), and this doublet is the lowest state also for 
the situation J± J||. Therefore the strong couphng 
limit of the asymmetric ladder with two spins s — 1, 
s = 1/2 is described as well by the triangular lattice 
depicted in Fig.Q in the same limit. The demand for 
in-triangle J|| to be large, in order to organize the effec- 
tive spin-1, is relaxed in this limit, and one can consider 
the situation with the uniform value of the exchange Jy 
along the whole leg. 

The presence of the anisotropy term —Dsfs2 in 129|) is 
a negligible effect in the described picture. Indeed, this 
term translates into a single-ion anisotropy of the triplet 
state, (5^)^ and the application of the formulas © 
shows that it is only the higher \Q) state, which becomes 
split accordingly, ~ D{S^)^. 

Let us discuss the analog of Eq. for the Hamil- 
tonian (|29|l . In the limit Jj_ 00 the effective inter- 
action is given by (jSJ with cr" acting within \Dl) and 
Jeff {00) = The analog of ^ for eq. ^ reads 



3 



(31) 



The second order in corresponds to transitions to 
higher \D0) and \Q) states. After some calculations we 
find 



J^ffiW ^ J^ffioo) 

+ei/37 -^ef fi°0)J]ff(.°0) 



(32) 



1 

3TI 



5 J_L + 2 J 



In H32|l the second term is negative for J± > J, in con- 
trast to 10. It means that the value of J^ff and the 



relative exchange anisotropy J^^ j / J; 
decreasing of Jj_. 

Thus we conclude that the strong coupling limit of 
the the trinagular ladder is described by the Hamiltonian 
© with Jeff given by Similarly to ((T^ . the in- 

leg correlations are slowly decaying ferromagnetic ones, 
although their modulation is different due to eq. H31|l 
instead of Q). 



B. Weak rung coupling, LSWT analysis 

For the case of ferromagnetic exchange with the easy- 
axis anisotropy we employ the spin- wave formalism. The 
situation is complicated by the existence of six spins in 
a unit cell. As a result, the quadratic Hamiltonian is 
represented as 6 x 6 matrix. The Hamiltonian for the 
interaction along the leg is standard, while the rung ex- 
change needs some care. 

Consider first two quantities Ai and Bi referring to ith 
site on the upper and lower leg, respectively. If Ai and 
Bi are coupled by the triangular rung exchange, Fig. |2 
then we have an expression 

J± 51 (^3J + ^3j + l)B3j + A3j+2{B3j + i + B3J+2), (33) 



with one term in the sum H33|) describing the coupling 
in the unit cell, and three-site periodicity of the overall 
structure. Going to Fourier components we have 

J2 + B-,+Jg + S-^-k/:,), (34) 



with K — 27r/3 and 

9,^lj±{l + e^') /, = -ij^(e- + e^?— ). (35) 

Particularly, for Aj = A = est, we obtain 

y^AiiBo + B^ + B^,). 

These preliminary notes show that the rung interaction 
hybridizes the magnons with the wave vectors q and q±K. 

The LSWT Hamiltonian H = J^k ^I-fffc*fc is obtained 
as a matrix defined for the vector 

The Green function Gij{k,t) — 
-i0{t){['i'i{k,t),^]{k,O)]) takes the form 
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where ojq = J| — J\\ cos q is the magnon spectrum for 
isolated chains, easy-axis anisotropy J| — J\\ > is as- 
sumed. The new spectrum is determined from the equa- 
tion (iei[G(g, — 0. The last equation amounts to 
the third-order polynomial in w^, which can be subse- 
quently solved. 

In order to analyze the lowest energies in the spec- 
trum, w ~ at g ~ it is sufficient to deal with smaller 
matrices. It can be shown that in this case one may con- 
sider almost degenerate 2x2 block formed by second and 
fifth lines (columns). The asymptotic expressions for the 
energies obtained this way coincide with those obtained 
directly from 

This simplified analysis can be also performed for other 
cases of the in-chain exchange anisotropy and rung ex- 
change. Particularly it is useful when the analytic treat- 
ment of the spectrum becomes problematic. For instance, 
the full LSWT consideration of the easy-plane D <Q for 
(I29|l amounts to the analysis of 12 x 12 matrix Green's 
function, while the simplified treatment reduces the cal- 
culation to bi-quadratic equation. 

In the subsequent equations of this Section, the rung 
exchange Jj_ appears with a prefactor 2gq=Q = 8/3. This 
prefactor is conveniently incorporated into the quantity 

Ji = 8Jx/3, (37) 

which is used below. Thus we interchangeably call Ji 
and J± as the rung coupling value. 

For the small anisotropy D, and |Ji| ^ J we find the 
following asymptotic expressions. 

i) Easy-axis, D > 0, AF sign J± > 0. Doubly degenerate 
gapful mode. 

ei,2,k ^ IVi'^D + Jk^){2D + 2Ji + Jfc2). (38) 

ii) Easy-plane, Z? < 0, AF sign J±^ > 0. One gapless, one 
gapful mode. 



ei,fc - -VJk^i^\D\ + 2Ji + Jk^), 

e2M ^ iV(2|^| + JA:2)(2Ji + JP). (39) 

iu) Easy-axis, D > 0, FM sign J± < 0. Two gapful 
modes. 

ei,fe~i? + ijfc^ s2,k^D+\Ji\ + ^Jk^. (40) 



Easy-plane ^ Easy-axis 







J± 




One gapless, 


one gapful 




Gapful, doubly degenerate 


Ai = 0, A2 


= 




Ai = A2 = V-Dp + Ji) 








D 


One gapless, 


one gapful 




Two gapful modes 


Ai = 0, A2 




-\Ji\) 


Ai=£l, A2 = £l + |Ji 



FIG. 3: The character of dispersion in different domain of 
parameters, the uniaxial anisotropy D and the rung coupling 
J±. 

iv) Easy-plane, D < 0, FM sign Jj^ < 0. One gapless, 
one gapful mode. 



ei,fc ~ -^Jk^i2\D\ + Jk^), (41) 

e2,fc ^ i V(2|D| + 2|Ji| + Jfc2)(2| Ji| + Jfc2). 

These results are similar with eqs. (|16|l . (|18|l . (|19|l . 
()20|l and will be compared below with the treatment by 
bosonization. Particularly they show that the low-energy 
dynamics of the triangular ladder is similar to the mixed- 
spin ladder of Sec. ^1 not only in the strong rung coupling 
regime, but also for the weak rung coupling. In Fig.|31we 
depict the character of dispersion in different domains of 
the small parameters, D, J±, according to eqs. (|38l) - (|41|l . 

C. Weak rung coupling, bosonization 

As discussed above, the spin-wave theory becomes in- 
applicable at larger J± when the role of quantum fluctu- 
ations grows. Instead, one may use the formalism which 
does not assume the average on-site magnetization and is 
suitable for spin-one-half chains. This formalism includes 
the Jordan- Wigncr transformation to spinless fcrmions, 
and the eventual continuum description with the use of 
fermion-boson duality described elsewherei^ 

This procedure is well defined for an easy-plane 
anisotropy, _D < in (|29|l , which is the case to be consid- 
ered in this subsection; the AF sign J± is implied. The 
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continuum representation of the spin operators in each 
of the chains reads as 



sf ^ e±'''(C + cos(7rx + 20)), 
~ TT^^dx'P + cos{'Kx + 2<l)), 



(42) 



with the omitted normahzation factors before cosines and 
constant C defined below. The bosonic fields (pi{r), 9i{x) 
are characterized by the chain index i and a continuous 
coordinate x. The Hamiltonian H = Hi + H2 + H± has 
a part for the noninteracting chains 



TTUK 



dx 

V 2 



-n 



2t:K 



(43) 



where i — 1,2 and 11.; = Tr^^dx0i canonically conjugated 
momentum to The form of H± in bosonization no- 
tation is discussed below. The general form of the leg 
Hamiltonian (|43|l prescribed by the bosonization proce- 
dure is complemented by the exact form of its coefficients, 
known from Bethe Ansatzii^ Denoting cos vrr/ = / J, we 
have 

l/K^2rj, u = J(sin7rr;)/(2 - 2?7). 
For \D\ J these formulas are simplified 



TTT] : 



m 
J 



K 



\D\J 



(44) 



and the coefficient C in becomes^^ ~ {nr])"^^ /8 ~ 
1/8. Note that the spinon velocity u in H44(l coincides 
with the one obtained by LSWT, eq. at Ji = 0. 

This unusual observation relates to the fact that the fer- 
romagnetic ground state is describable in classical terms 
of the total magnetization, and the semiclassical LSWT 
approach should work well in the nearly ferromagnetic 
situation even without LRO. 

The rung interaction H± couples different terms of the 
spin densities. However the interaction of the AF com- 
ponents, ~ cos(7ra;-|-20) in (|42|l is irrelevant in the renor- 
malization group sense ; moreover, the structure of H±, 
ea. (|34|l . shows that this interaction is absent in the lowest 
order. 

It is convenient to introduce the symmetrized combi- 
nations 0± = (0i±02)/V2, 8± = (6'i±^2)/V2. In terms 
of these, the relevant and marginally relevant terms of 
H± are 

H± ~ Ji J dx[C^ cos{V2e-) + 8-^/^0^0x9+ sm{V2e-) 

-f (2^)-2((a,^+)2 - (a,</>_)2)]. (45) 

The second term in ()45|l conies from the gradient expan- 
sion of (|34|l . Its inclusion however does not change the 
Gaussian character of the action for the field 6+ (see be- 
low). Integration over this latter field produces a contri- 
bution ~ {Jf / J) cos2^/26^ , which i) is less relevant and 



ii) has a smaller prefactor than the first term in H45|l . 
That is why we omit the term dx0+ sin (a/26'_) below. 

The remaining terms are combined into the Hamilto- 
nian of the form H = H> + H_ with 



dxi^^^n' 



= dx{ 



2 --+ 



2ttK^ 



-(9.0+)2), (46) 



2 - 2ttK^ 
JiC^ cosy26l_). 



(47) 



where 



U+K+ = u^K^ = uK ~ 7rJ/4, 
u±lK± = u/K±Ji/2t:-^{2\D\±Ji/2)/it. (48) 

Eqs. (|46|l . H48|) show that the mode 0+ remains gapless, 
and its velocity is increased with Jj_ , hence it corresponds 
to El mode in (jSlJl, ((T^ . 

The situation with the 0_ mode is more complicated. 
Two features are noted here, the instability of the Gaus- 
sian action at Ji > 4|Z?| and the appearance of the gap 
in the spectrum. 

Indeed, the scaling dimension of the operator cos \/2^- 
in (|47|l is l/(2if_) ^ 1 and the dynamics of 9_ mode 
is desribed by the sine-Gordon model in the quasiclassi- 
cal limit, with a large number of quantum bound states, 
or " breathers" . The gap A in the spectrum of 9- field 
is given by the mass of the lightest breather, which is 
roughly found by expanding the cosine term and rescal- 
ing the field 9- Q-j ^JkZ 



A2~2™_JiCVi^-- ^ 



1^1 -T 



(49) 



In the leading order in Ji this expresssion corresponds 
to the mode £2 in (|^ . ((T^ . The refined value of the 
gap can be obtained after usual scaling argument»i2*i^ or 
directly from the exact formulas ini&. The identification 
of our model parameters with those of Lukyanov and 
Zamolodchikov reads as 

[I = u_ J1CV2, = i'iK^)-^ 

and u_ stands for the overall energy scale. 
The gap (m in notation of^S,) is then found as 



a2 J 



and the spectrum becomes 
1 



4,. 



J 



1/(1-1/(4K-)) 



(50) 



-[\D\ + ^]Jk 

1 
2 



(51) 



The mean value of the cosine term is given by the ex- 
pression 



(cos 



I 16 J 



(52) 
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According to H48() . H52|l . the increase of Ji leads to the 
instability of the Gaussian action, which happens simul- 
taneously with the saturation of the quantity (cos •\/26'_) . 
Note that the similar situation was observed inpiS for the 
simple two- leg ferromagnetic ladder. 

For one chain, the breakdown of the Gaussian action 
happens at _D > 0, and corresponds to the transition to 
the ferromagnetic ground state. The average value of 
spin in this case becomes (s^) = n~^{dx(f>) = ±1/2. 

For a ladder, the discussed instability and the satura- 
tion of cosine term correspond to the saturation of the 
scalar product of spins in different chains, (sijS2j) (see 
below). It means that the spins in adjacent chains form a 
singlet state. The peculiarity of this phenomenon is the 
energy scale when it happens, J± \D\, following from 
the bosonization (cf>i^). 

This small energy scale is unusual and may be com- 
pared to the LSWT treatment. Successful enough for 
isolated chains, LSWT is in qualitative agreement with 
bosonization, regarding the increase of the spinon veloc- 
ity with J± for the symmetric mode 0+ as well as the 
gap value for the 0_ mode. At the same time, LSWT 
predicts the increase in the velocity of with J± at the 
energies higher than the gap value. The bosonization 
says the opposite. 

At this moment it is also instructive to consider the FM 
rung coupling Ji < 0. In this case the LSWT formulas 
(Pn|l . l|lT|l show again one gapless and one gapful mode, 
with the unchanged and increased velocities, respectively. 
The bosonization, H51|) . provides the similar picture, but 
says again about the collapse of the gapless 4>+ mode at 
\Ji\-\D\. 

It is worth to note here that the average cosine term 
(|52|l and the correlation length, eq. (|55|) below, does not 
show any peculiarities at Ji ~ \D\. 

A possible explanation for the above discrepancy stems 
from the observation that £2 mode (|39|l at Ji = \D\ at- 
tains the form 62 — \D\ + Jk^ /2. The region of linear 
dispersion of bosons, a cornerstone of conformal treat- 
ment, is lost here, which may be reflected by vanishing 
velocity in the bosonization treatment. Note also, that 
eq. I|39|) shows the roughly linear gapful spectrum upon 
the further increase of the rung exchange, Ji \D\. 
This feature should assumably be valid in the corrected 
bosonization treatment. 

We suggest here that the action l|T7|l should be comple- 
mented by the irrelevant terms, usually dropped in the 
infrared limit. They come from the consideration of the 
lattice Hamiltonian and are of the structure 



(53) 



with a the lattice spacing. The appearance of these terms 
is most easily observed by the consideration of one-chain 
XY model. In terms of the Jordan- Wigner fermions ip , 
one has the tight-binding fermionic spectrum, cos q. Near 
the Fermi points q — ±tt/2 the leading terms in the ex- 
pansion of the fermionic dispersion are the linear and 
cubic terms. The linear-in-g fermionic term,~ %p^dxfp, 



transforms into {dx'p)'^ in the bosonic language, and the 
cubic term ~ ijj'^d^ip attains the form H53|l. Omitting the 
unknown coefficients ~ 1 and denoting dxip- = 4>'^ 
etc., the new Hamiltonian H47|l is then schematically writ- 
ten as 



J 



2 - 



\D\ - Ji 



4,' ^2 



J, 



J, 



r+^('^''r+T('^-r (54) 



Let us consider first the case Ji = 0. The interac- 
tion term {dx4>)'^ may be discarded in the infrared action, 
and the quadratic term (9^(^)^ modifies the spectrum, 
e| ~ \D\J}t^ + J'^k^, so the spectrum may be regarded 
as linear only at fc < y/\D\/J. The latter estimate is 
in accordance with the LSWT formulas dl). The 
dynamical correlation function becomes 



{dx4>-,dx4>-)ku 



Jk^ 



;2 - J|D|fc2 _ J2fc4 



which leads to the estimate for the average square 
(correlation function a,t x ^ t — 0), {{dx<j)-Y) ~ 
/jj dkk/ yjk"^ + \D\/ J. The latter quantity is defined by 
large fc 1 and shows that the fluctuations are strong, 
{{dxii-Y) ~ Ij should be expected from H42|) for the 
fluctuating spins without LRO. 

Consider now the case Ji 7^ in (|54|) . In the regime 
with the negative coefficient before {dx(f>-)'^, the inter- 
action term (9^0)'* stabilizes the action against the di- 
vergent static mean value 3x4'- ■ The usual recipe here 
is first to determine the variational static solution to 
the above Hamiltonian letting n_ — — 9^, see, 
e.giiS and references therein. The trivial classical so- 
lution is the doubly degenerate vacuum dxip^^"^ = po ^ 
±1/ ( Ji — \D\)/ J. The spectrum of fiuctuations around it 
is well-defined with the velocity ^ J( Ji — \D\) and the 
Luttinger exponent ~ J/{Ji ~ The short-range 
fluctuations are still determined by the quadratic part of 
the spectrum, Sk ^ Jk^, and the average square of the 
fluctuations is similarly estimated, {{dx4>- — dx4>-^)'^) ~ 
1 . It shows that the amplitude of the quantum fiuctua- 
tions exceeds the distance between the vacua, 2po, which 
makes the choice of the classical vacuum dubious. 

The refined analysis reveals the existence of multi- 
soliton classical solutions to H54|l . Variating the static 
Lagrangian over (/)'{x) and letting 4>'{x) = pofiu) with 
y — pqx we obtain an equation —d^f/{dy)^ = f — f^, 
which allows a solution of the form / = Q;isn(Q;2y, k) 
with sn(?/, k) the Jacobi elliptic function and k the elliptic 
index. It leads to the iV -soliton solut ion for classical vac- 
uum, cJaj^^"* — ±po-\/2k2/ (1 + k'^)sii{xpq/\^1 -t- k^, k), 
with the soliton density N/L = po/(2Vl + ^^^(k))^^ 
In the limiting case k = 1 one has one soliton, dx<p'^^'^ ~ 
po tanh(a;po/V2)- The difference in the classical energy 
between these A^-soliton solutions and the above trivial 
vacua is estimated as ~ NJp^, i.e. a small quantity 
a,t N ^ 1 <^ PqL, as compared to the classical energy 
^ LJpq. The full analysis of the problem should hence 
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include the summation over the A''-sohton solutions. The 
existence of the quantum gap, ~ Ji{Ji — \D\), ex- 
pected from the Jid^ term in H54|l . only adds to the 
complexity of this problem, which should be dicussed 
elsewhere. One can only observe here that the neces- 
sity of summation over the classical vacua provides the 
absence of the staggered magnetization along the z— axis, 
associated with the non-zero classical dx(t>-- 

Knowing the spectrum and the Luttinger exponents, 
one can use the principal advantage of the bosonization 
in evaluation of the correlation functions. These correla- 
tions are discussed in the next section upon the assump- 
tion of the weaker coupling, Ji < \D\. 

D. Correlation functions 

The spectrum of H47|l consists of one gapless and one 
gapful mode. The gap A corresponds to a finite correla- 
tion length 

^ = u/A^ (55) 

separating domains of different behavior of the correla- 
tion functions. The transverse spin correlations in one 
chain, j = 1,2, have the formSP 

(4oS,>) - r-i/4^^e(^«('-/«)-^«(«/«))/(4^-), (56) 

with Kq{x),Ki{x) modified Bessel functions and a the 
lattice spacing. At shorter distances, r < ^, this expres- 
sion becomes (s^s^) ^ j,-i/4A:+-i/4_fs'_ ^ while at large 
r > ^ one has (s^s^) ^ ^-i/4-ft:+^-i/4i<'_ ^ rpj^g interchain 
correlations are 

{s+QS2r) ^-l/4if+g-(ifo(r/0+^oW5))/(4if-)^(-57) 

with the behavior _r-i/4if++i/4X_^-i/2K_ 
_^-i/4K+^-i/4K- g^Qj-ter and larger distances, 
respectively. Hence the interchain correlations decay 
faster beyond the scale r ^ 

The longitudinal correlations are obtained in the form 

{4,o4,r) - ^ + ^K,ir/0, (58) 

(4o4.) - ^-^K,ir/0, (59) 

which shows particularly that at r < ^ the interchain 
correlations are of the AF character. 

The parameters K^^ ,u+, A,(^~^ increase with J^. 
Hence the transverse correlations decay faster at larger 
J±. We argued above that in the strong-coupling limit 
J_L ^ oo one deals approximately with the AF Heisen- 
berg chain situation, wherein (s^s^) ^ r^^. Compar- 
ing it with H56|) . (|57|l one may conclude that K^^ should 
reach the value 1/4 in the strong coupling regime J± ~ J. 
Actually it is not so simple, since the derivation of H47|) 
assumed > 1, and other terms of the rung interac- 
tion become important at smaller K^. As a result, one 



expects that the increase of J± eventually changes the 
structure of the effective low-energy action. 

Summarizing, we show that the "triangular" model 
of this section is equivalent to one of Sec. for the 
strong rung couplings. Further its dynamics is similar 
to one of the mixed spin ladder also for the weak rung 
coupling, as shown by LSWT approach complemented 
by the bosonization. Two latter techniques reveal cer- 
tain shortcomings in the decription of the situation, as 
LSWT becomes formally inapplicable without LRO and 
the bosonization becomes unstable at the level of the 
Gaussian action. 

Working in the close vicinity of the FM point in the 
parameter space {D,J±), we observe that the transi- 
tion to the FM ordered phase is of the first order at 
the line D 0, J± =0. At non-zero AF values of J±, 
this transition becomes the second-order one, at the line 
^ D > 0. Approaching this transtion line from above, 
J± > J^, one should observe the divergence of the cor- 
relation length and vanishing critical exponents of the 
correlation functions. 

Combining the results of Sec. Inland Sec. lIIII we conjec- 
ture that the crossover from the weak to strong rung cou- 
pling regime for the isotropic situation, = 0, is charac- 
terized by the absence of the long-range order and gapless 
character of dispersion. Increasing the AF value of J± on 
has £fc ~ \/irj±^k until J± < J and Sk ^ Jk at J± > J. 
This form of dispersion takes place at fc < ~ \/ J±/J- 
The correlation functions are of the form (sg s") ~ r~'^ 
beyond the correlation length, ^, with 7 ~ \/J±/J at 
J± ^ J and 7=1 otherwise. 

Note that this non-universal behavior of the critical 
exponent 7 characterizes the isotropic gapless situation 
and should be contrasted to the well-studied case of a 
gapless XXZ chain. ^ In the latter case one has different 
exponents for different spin projections a, with certain 
relations between them, e.g. 7^:72 = 1/4. 

IV. ORDER PARAMETERS 

A. String order parameter vs. scalar product 

In the pape»2i (see also^SiSi) a model of a symmetric 
AF Heisenberg ladder of spins s = 1/2 was considered. 
Particularly, the authors discussed the string order pa- 
rameter (OP), which was associated with the topological 
OP introduced earlier— for the spin-1 chain. In fact, the 
discussion by Shelton et al?^ for non-zero AF rung ex- 
change Jj^ can be reduced to the observation that the 
scalar product S1S2 on the rung assumes the non-zero 
value. 

Let us characterize each state of two spins on a rung j 
in terms of singlet \Sj) and triplet \Tj). The ground state 
\G) of the whole ladder has a component comprised of 
all rung singlets, \Stot) = ®j\Sj). It is clear that for the 
case of extremely large AF rung exchange the weight W 
of \Stot) in \G) is unity. One expects that for moderate 
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AF J|| this weight W is finite. Consider now the 

spin product on jth rung — 4s"]^s"2 — expi7r(s"]^ + 5^2) 
with a = x,y, z, which may be represented as 

-Ki'h = -Ps,, + (1 - T'^.^Oe"^" (60) 

with Vsj projecting onto the jth singlet and S°' spin-1 
operator for the jth triplet. Note that the presence of 
Vsj makes (|6U|) different from the operator e'^'^J used 
by den Nijs and Romnielse in their discussion^** of the 
spin-1 chain. 

Indeed, the "string" operator ]Xj=i{~^^'j,i^'j,2) its 
ground-state expectation value contributed by the weight 
of the \Stot) state. This partial contribution is equal to W 
and does not depend on the distance [n — l). Particularly, 
the expectation value of the scalar product —Asj iS.j 2 = 
iVsj — 1 has a contribution 3W from \Stot) state. In 
bosonization notation we have, 

(sj_iSj_2) ^ (cos\/20_) — (cosa/8(/)+) + (cos \/8(?!)-). 

For the AF signs of J^\,J±, considered inSi, first two 
cosines in the latter expression have nonzero values. 
Some inspection shows that these values correspond to 
ones reported inSi for the infinitely long string OP. 

Hence we conclude that the string OP discussed in^^'^'^ 
for AF rung interaction can be identified with the scalar 
product of spins on a rung and measures the weight W 
of the total singlet in the ground state. It should be 
stressed, that our above arguments are not applicable 
for the FM rung interaction, when the lowest rung state 
is triplet. In this latter case the string order parame- 
ter discussed ii^SLSSiS^ is the appropriate description and 
cannot be reduced to a local scalar product. 

Clearly, the non-zero average value of the scalar prod- 
uct is disconnected from the appearance of the on-site 
magnetization, as discussed below. 



B. Asymmetric ladders 

Applying the same type of consideration to our above 
systems, we can say, e.g., that for the mixed spin ladder 
the order parameter is the average value of the scalar 
product on the rung pj — SjSj. It assumes two values, 
Pj = —1 and 1/2 for the rung doublet and quadruplet, 
respectively. For J± = all these six states have the 
same weight, resulting in (pj) = 0. With the increase of 
AF rung exchange, (pj) saturates into —1 value. 

Similarly, for the triangular ladder, one considers the 
combined scalar product pA = (si -I- 82)33, see Eq.ljJOJ). 
This quantity takes three possible values pA = —1, 0, 1/2 
in the states \D1), \D0), \Q), respectively. Increasing J±, 
the |D1) state becomes favorable, with {pa) — 1- 

Notice, that the discussed order parameter is bilin- 
ear in spins, independent of the in-leg spin exchange 
anisotropy and does not imply the ordering of individ- 
ual spins. 



The spin ordering in a proper sense depends on the 
sign of the uniaxial anisotropy. Particularly, in the case 
of the easy-plane anisotropy, both the weak and strong 
rung coupling regimes correspond to XXZ model in the 
absence of LRO. Therefore one does not expect the spin 
ordering here. 

The case of the easy-axis anisotropy can be analyzed 
for the mixed spin model. We showed in Sec. that 
the LSWT, applicable for isolated chains, fails for the 
intermediate J±. At the same time, the strong cou- 
pling Hamiltonian ||SJ| is the AF easy-axis XXZ model. 
It means the appearance of non-zero staggered magne- 
tization for the effective spins (t| in Scaling esti- 
mates (see e.gM.) show that (cr|) - {-iy{D/J)°' with 

a — (t^I'^^ J jD. This exponentially small value of the 
order parameter for the effective Hamiltonian (jSJl trans- 
lates into the corresponding values for initial spins ac- 
cording to l@J. Note that the average spins in one leg 
are aligned in one direction, but due to the difference 
in their contribution to the rung doublet state, eqs. I0J, 
(|31|l . both the uniform and the staggered magnetization 
is present in the system^ 

V. CONCLUSIONS 

We demonstrate above that the mixed spin ladders and 
triangular ladders with the ferromagnetic coupling along 
the legs are generic models for description of a transi- 
tion from the classical (ferrimagnetic) to quantum (an- 
tiferromagnetic) regime. The individual legs with the 
isotropic Heisenberg exchange show the classical ground- 
state and their dynamics is well described by the quasi- 
classical spin-wave theory. Turning on the AF rung cou- 
pling introduces strong fiuctuations, which destroy the 
long-range order and eventually make the system equiv- 
alent to the quantum AF spin s = 1/2 Heisenberg model. 

We showed that in a large domain of parameters for 
these ladders the spin wave theory, although missing cer- 
tain features caused by quantum fluctuations in one di- 
mension, is still quite instructive for the qualitative deter- 
mination of the spectra, which allows the further compar- 
ison with more sophisticated methods. The refined anal- 
ysis of the spectrum and correlations by the bosonization 
technique complements the investigation of the "quan- 
tum" regions of the phase diagram. As a result, the 
unified description of the model becomes possible, partly 
including the complicated crossover region from the weak 
to strong rung coupling limit. 

We argue that for the isotropic spin exchange this 
crossover is characterized by the gapless spectrum with 
(spinon) velocity ~ \f77Z. The vanishing velocity at 
J_L = corresponds to the first-order phase transition to 
the ferromagnetic state. The asymptotic decay of corre- 
lation functions is described by unique critical exponent, 
7 ^ yJ^JJ ^ for £^11 three projections of spin. This type 
of behavior makes the mixed spin ladder in its crossover 
regime quite distinct from the AF s = 1/2 Heisenberg 
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model, which would be very interesting to verify by in- 
dependent, e.g. numerical, methods. 
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